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A simple model of species origin resulted from dynamic features of a population, solely, is devel- 
oped. The model is based on the evolution optimality in space distribution, and the selection is 
gone over the mobility. Some biological issues are discussed. 
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I. INTRODUCTION 

Modelling of the dynamics of biological communities 
is essential for further understanding both mathematics 
[l|,[3,f3, 4], and biology [i,[i,0,ll- A lot have been done 
in mathematical ecology and mathematical population 
biology since the pioneering works by V. Volterra [9] and 
A. Lotka [13] ■ Yet, there is no comprehensive, steady 
and solid general theory of the dynamics of biological 
communities. In spite of the implementation of the most 
general and outstanding theorems of the natural selection 
[111 . [T3 , quite a number of problems still await a 
researcher to deal with. 

Modelling of evolution processes seems to be the only 
way to figure out various important and significant de- 
tails in evolution theory. Since Darwin's famous work 
(T^ . evolution theory evolved heavily, itself. A lot 
have been done in order to introduce the mathemati- 
cally based methodology into the selection theory; Hal- 
danc'outstanding pa pers made the most successful start- 
up, in this direction [l5l. [l^. 

Here we propose a simple model of a species origin 
resulted from the dynamics of a population, solely. An 
origin here means a dissociation of an originally uniform 
population into two subpopulation distinctively differing 
in the mobility. This discretion is understood as an ap- 
pearance of a polymodal distribution of the beings over 
the mobility character. To begin with, we should con- 
sider, in detail, the model of optimal migration. 



II. BASIC MODEL OF OPTIMAL MIGRATION 

Modelling of spatially distributed communities and 
populations is quite a problem. Basically, the approach 
based on reaction -i- difFusion methodology is used, for 
modelling of those communities. The diffusion approach, 
being quite attractive from mathematical and physical 
point of view, has nothing to do with any real system (or 
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FIG. 1: A distribution resulted from the selective pressure. 
Here is the distribution in a single site shown; there are no 
beings in other one. 



even any system just pretending to be relevant to a real 
one). That is the diffusion, that makes the problem here. 

Diffusion approach gets very strict and absolutely un- 
feasible constraints on individuals under consideration: 
they must transfer in space in random and aimless man- 
ner. There is no one species going this way; even bacteria 
in continuous cultivation systems control, to some extent, 
their location. The hypothesis underlying the diffusion 
methodology could be improved in neither way. 

Thus, the approach based on the principles of evo- 
lutionary optimality [111, [H, [13 opposes the diffusion 
methodology; originally, these principles have been in- 
troduced by J.B.S. Haldane. From the point of view of 
the dynamics of spatially distributed communities, these 
principles force the individuals to migrate in the manner 
to increase the net reproduction (NRG). NRG is an 
average number of descendants (per capita) determined 
over a sufficiently long generations chain. In brief, a 
model based on the the principles should look like the 
following. Gonsider a population occupying two sites; a 
transfer from site to site makes a migration. No other 
movements, inevitable in any real situation, would be 
taken into account; moreover, we shall suppose that such 
movements have no impact on the dynamics. 

Suppose, then, that the dynamics of each subpopula- 
tion follows the Verchulst's equation: 
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FIG. 2: The case of suppressive selection: the beings of the 
pecuhar mobihty class {p = 0.8) are eliminated. 



and 



Mt+i = Mt ■ (c - dMt) 



(lb) 



Here Nt {Mt, respectively) is the subpopulation abun- 
dance determined at the time moment t. We shall con- 
sider the dynamics in discrete time; continuous time con- 
sideration is possible, as well, while it bring nothing new, 
but the severe technical difficulties. 

We have chosen the Verchulst's equation due to its 
universality The linear functions in ^ are NRC, in 
relevant sites. Thus, a migration from one site to another 
starts up, as soon, as the living conditions (measured as 
a part of NRC) becomes worse here, in comparison to 
similar ones there, with respect to the transfer cost p, 
< p < 1. This latter may be considered as a probability 
of a successful transfer from a station to another one, 
with no damage of a further reproduction. Evidently, an 
additive pattern of a transfer cost is the simplest one: 

P = Pout + Pin + d , 

where Pout is the cost of a successful leaving of the site; 
Pin is the cost of a successful intrusion into another site, 
and d is the pure transfer cost. 

We shall suppose that the individuals under consid- 
eration are globally informed; in such capacity, it means 
here, that all the parameters (including the transfer cost), 
as well, as the abundances in each site are known to them. 
The 
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migration takes place, as 
b- N <p- [c- d- M) , 
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FIG. 3: An occurrence of several subspecies differing in mobil- 
ity. It is clearly evident, that some classes of p are eliminated 
Upper figure corresponds to 
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site, and lower one to 
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site. 



and the reverse migration takes place, when 
c~d- M <p - {a~b-N) . 



(3) 



The number of migrants must equalize the inequality ^ 
(or the inequality ([3]), respectively): 

a-b-{N-A)=p-{c-d-{M+p-A)) (4) 



c-d- {M - A) ^p- {a-b- (N + p- A)) , (5) 

respectively. Since the number of emigrants may not ex- 
ceed the total abundance observed within a station, then, 
finally, the migration flux A is determined according to 



A = min i 



or 



A = min <^ M, 



pc~ a + dNt — pdMt 
b+p^d 

pa- c + dMt - pbNt 



d + p'^b 



(6a) 



(6b) 



Here (|6ap corresponds to 
versa. 
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migration, and vice 



III. MODIFICATION OF THE BASIC MODEL 

Suppose, a population (in both sites) is divided into 
K classes, with respect to the mobility character. This 
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FIG. 4: Another example of disruptive selection. 

latter is measured as transfer cost peculiar for each class: 

< pi < P2 < ■ • • < PK-i < < 1 ■ 

Let, then, A^^-''' {M^^\ respectively) be the abundance 
of of a subpopulation at j-th class. The dynamics of a 
subclass (at a subpopulation) is determined according to 
the equation: 

iV(^,=7Vj^").|a-6.l|j7Vi^)| ; (7a) 

M« =M(^").|c-d.lf:M«| . (7b) 

Migration rule keeps the same, as ^ and ([3]), for each 
class, with pj peculiar for the class. So, the migration 
fluxes are determined independently, at each mobility 
class (i. e., for eachpj), the beings relocate themselves be- 
tween the sites, and finally the reproduction takes place. 

IV. RESULTS AND DISCUSSION 

We simulated the dynamics of the system ([D - El) , for 
various parameters. Obviously, the system is symmetric: 
nothing will change, if one simulate the dynamics origi- 
nally for a — a,b = P^c = "f and d = 6, and then do it for 
a — "f,b — 9,c = a and d = p. This symmetry allows to 
decrease the dimension of the system. Three parameters 
control the system: these are a, c and A = b/d. 

We have simulated the dynamics of the system (IT]-[7]) 
for A = 1, A = 1/2, A = 1/4 and A = 1/10. The simu- 
lation was run in the manner ro scan the parameters a 



and c area searching the stable polymodal distribution 
of the the beings over the transfer cost values. More- 
over, in case of A = 1, there is a symmetry against the 
permutation of the parameters of a and c. 

Further, we show several figures illustrating the ob- 
served distributions. Fig.[T] shows the distribution ob- 
served at A = 1/4 and a — A, c = 0.1. The observed 
distribution is a typical pattern, for the case of so called 
selective pressure: the abundance of the classes with de- 
creasing mobility grows up monotonously. The figure 
shows quite rare situation, when the beings with low mo- 
bility (small p) take an advantage under the selection 
pressure. 

Figure [5] shows an interesting case of the suppressive 
selection: a peculiar mobility class {p = 0.8). This kind 
of the distribution might be considered as a partial case of 
disruptive selection. Similar situation is shown in Fig.|4l 
again only one site is shown, since the other one is not 
inhabited. This figure is obtained for A = 1/2, and a = 
0.1, c = 8.7. 

Figure [3] shows the situation that is clearly understood 
as the dissolution of originally homogeneous (and con- 
tinuously distributed over the mobility factor p) popula- 
tion into two discrete subspecies. Some classes are com- 
pletely eliminated due to the selection. Here only few 
classes have survived; the survival of the most mobile 
classes seems rather natural, while the pattern shown in 
Fig. [1] puts on a hypothesis that there are some param- 
eters vales that would yield the survival of the slower 
beings, while the faster ones would be eliminated. Fi- 
nally, Figure [3] shows the situation of the elimination of 
the peculiar classes, while the majority of them survive. 

The paper aims to illustrate the occurrence of the se- 
lection based on the mobility of organisms, in case of a 
non-random migration. Actually, such selection is not a 
trick: according to the fundamental theorems on the se- 
lection, everything that is inherited, must be processed 
through the selection [ll|, [l^, [H, [l3| ■ This fact explains 
the results observed at the simulation experiments. 

Obviously, one can argue, that there are no species that 
are present by the beings with all possible diversity of a 
parameter (the mobility p, in our case). That is right; 
here we just went the way very peculiar for a modelling 
in biology. We have changed the potential diversity of 
the beings with mobility p that may appear within a 
population due to mutation for the actual diversity, as if 
all the mutations already have been done. 
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